Electromagnetic corrections in 77 ^ Stt decays 
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Abstract. We re-evaluate the electromagnetic corrections to Sn decays at next-to-leading order in the 
chiral expansion, arguing that effects of order e^(m„ — nid) disregarded so far are not negligible compared 
to other contributions of order times a light-quark mass. Despite the appearance of the Coulomb pole 



m T] 



TT^TT tt" and cusps in 77 ^ Sty" , the overall corrections remain small. 
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1 Introduction 

The decay 77 ^ Stt is particularly interesting because it 
is forbidden by isospin symmetry. While the 77 meson has 
isospin / = 0, three pions with zero angular momentum 
can only couple to / = 1, and thus the decay can only 
happen via isospin breaking AI — 1 operators. In the 
Standard Model, there are two such sources of isospin vi- 
olation, on the one hand strong interactions from 



rrid - niu 



[dd — uu){x) , 



(1.1) 



which are proportional to the light-quark mass difference 
vrid — ruu, and on the other hand electromagnetic interac- 
tions that are proportional to the electric charge squared 
from 



dyD^^''{x-y)T{3^{x)Uy)) , (1-2) 



where D'^''{x—y) is the photon propagator and j^(x) is the 
current density containing the charged fields of the theory. 
Sutherland and Bell showed by using soft-pion techniques 
that the electromagnetic contribution at tree level is much 
too small to account for the observed decay rate [1,2]. The 
rj decay is therefore very sensitive to nid — rriu and hence 
it potentially yields a particularly clean access to the de- 
termination of quark mass ratios. The strong tree-level 
amplitude was subsequently studied using current algebra 
and partially conserved axial-vector current (PCAC) tech- 
niques [3, 4], but the decay width turned out to be off from 
the experimental value by a factor of a few. Later PCAC 
and current algebra were generalized and cast in a mod- 
ern form in the framework of chiral perturbation theory 
(ChPT) [5-7]. Gasser and Leutwyler (GL) subsequently 



calculated the strong contribution at one-loop level [8]. 
They observed large unitarity corrections due to strong 
final-state interactions, but their value still differs from 
experiment by a factor of 2. Corrections beyond one loop 
were studied using dispersive methods [9,10], unitarized 
ChPT [11,12], and finally with a complete two-loop cal- 
culation [13]. All of these find considerable enhancement 
of the decay widths compared to the one-loop calculation. 

Baur, Kambor, and Wyler (BKW) [14] studied correc- 
tions to Sutherland's theorem by evaluating the electro- 
magnetic contributions in 77 ^ Stt at one-loop level, using 
an extension of ChPT including virtual-photon effects [15], 
but they found them to be very small. The motivation 
for reconsidering these electromagnetic corrections at one- 
loop chiral order in the present work hinges on the fact 



that Ref. [14] neglects terms proportional to e^iirid - 



(such terms have recently been considered for r/ — > 37r*' m 
Ref. [16]), arguing that these are of second order in isospin 
breaking and therefore expected to be suppressed even fur- 
ther. However, by restricting oneself to terms of the form 
e^m (where m = [rriu + md)/2) and e'^nis, one excludes 
some of the most obvious electromagnetic effects: real- and 
virtual-photon contributions, as well as effects due to the 
charged-to-neutral pion mass difference (which is predom- 
inantly of electromagnetic origin), both of which scale as 



'(md- 



These mechanisms fundamentally affect the 



analytic structure of the amplitudes in question: in the 
charged decay channel rj 7r+7r^7r°, there is a Coulomb 
pole at the boundary of the physical region (at the tt+tt" 
threshold), while in the neutral decay channel rj — > Stt^, 
the pion mass difference induces a cusp behavior at the 
7r+7r~ thresholds (compare e.g. Refs. [17,18]). This cusp 
encodes information on tttt scattering in principle in much 
the same way as the decay — + tt'-'tt'^tt^ , which has been 
established as a new means for a precision determination 
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of the scattering lengths combination gq — a2 [19-25]. In 
contrast, the corrections identified in Rcf. [14] arc all poly- 
nomials (due to counterterms) or quasi-polynomials (due 
to kaon loop effects) inside the physical region. Further- 
more, a soft-pion theorem [1] guarantees that at the soft- 
pion point, these corrections are of order e^m only (and 
not O(e'^ms)); hence the relative suppression of the ne- 
glected terms is of the order of (mrf—m„)/(mrf-fTO„) ~ 1/3 
and therefore not a priori small. 

This work is organized as follows. In Sect. 2 we give a 
short introduction into the formalism of chiral perturba- 
tion theory with virtual photons and mention some special 
topics that are relevant for our investigation. The decay 
amplitudes for both channels are derived and displayed in 
Sect. 3, and in Sect. 4 our numerical results are presented. 
Finally, we conclude with a summary and an outlook. Var- 
ious technical issues are relegated to the Appendices. 



2 Formalism 

A particularly successful approach to describe the interac- 
tions of hadrons at low energies is to construct an effective 
field theory which encodes the infrared behavior of QCD. 
Based upon the approximate chiral symmetry of Cqcu one 
can identify the lightest particles in the spectrum with the 
pseudo-Goldstone bosons induced by spontaneous break- 
ing of this symmetry, and their interactions can be writ- 
ten as an expansion in small momenta and light-quark 
masses. This effective theory of the strong interactions is 
called chiral perturbation theory (ChPT) [5 7] . It can be 
generalized to include electromagnetic effects systemati- 
cally [15] (see however Rcf. [26]). There are several good 
introductions to ChPT; see for example Refs. [27-29]. 



2.1 Lagrangians 

The Goldstone boson fields are collected in the field U 
according to 



U = exp 



V2 



JlA. _|_ 


V8 

ye 


7r+ 




K+ 


TT 




•"■3 1 


V6 


K° 


K- 




^0 


V6 



(2.1) 

The effective Lagrangian for the mesonic sector invariant 
under chiral symmetry reads at leading order (LO) 0{p^) 



^(2) _ ^0' 



(D'^U^D^U + x^U + U^X) , (2.2) 



where (...) denotes the trace in flavor space. The external 
(pseudo-) scalar sources are conventionally combined into 



X = 2Bo{s + ip), s = M + ... 



(2.3) 



incorporating the quark mass matrix A4 = 
diag(m„, TOd, TOs). Fq is the pion decay constant in 
the chiral limit, and Bq is related to the chiral quark 
condensate. 



The terms of the Lagrangian at next-to-leading order 
(NLO) 0{p*) needed in the following are given by 

+ U{D''U^D^U){x^U + U^x) 
+ L5{D'^U^D^Uix^U + U^x)) 
+ Leix^U + U^x)^ + Lrix^U - U^x? 
+ Ls{x^Ux^U+U^xU^x) ■ (2.4) 

(2) 

Loops with vertices from the LO effective Lagrangian Cl^^ 
generate infinities that can be absorbed by a renormaliza- 
tion of the introduced NLO low-energy constants (LECs). 
For this purpose one defines renormalized low-energy con- 
stants 

Li = riX + Ll{fi), (2.5) 

where A contains a pole in d = 4 space-time dimensions, 
and the coefficients 1^ are given in Ref. [7]. The renor- 
malized coefficients depend on the scale /U introduced by 
dimensional rcgularization. All physical observablcs are 
finite and scale independent, which serves as a check in 
explicit calculations. 

Following Ref. [15] to include electromagnetism in the 
framework of ChPT by adding virtual photons as addi- 
tional dynamical degrees of freedom, the full NLO effec- 
tive Lagrangian can be written as 

£eff = C^l + 4S + -C^S + 4S + . (2.6) 

The new local electromagnetic interactions in the La- 
grangians cim, jC-im contain the quark charge matrix 
Q = ediag(2, — 1, — 1)/3. At leading order the local in- 
teractions stem from the following additional effective La- 
grangian 



4'i = -V''^^. + C(QC/QC/t), 



(2.7) 



where F^^ is the electromagnetic field strength tensor, 

and we have omitted the gauge fixing term. The new low- 
energy constant C determines the purely electromagnetic 
part of the masses of the charged mesons in the chiral 
limit. For convenience one often defines the dimensionless 
constant 

Z=^, (2.8) 

which we will use in the following. 

At next-to-leading order the terms of the additional 
Lagrangian relevant in the following read [15] 

4^ = F^{{D^'U''D^U) [Jri(Q2) + /^2(Qf/Qf/t)] 

+ Ks{{D^^U^QU){D^U^QU) 

+ {D^UQU^){D^UQU^)) 
+ Ki{D''U^QU){D^,UQU'^) 

+ Kq{DI'U'^D^UQU'^QU + Df^UD^U^'QUQU^) 
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+ K^iix^U + U^x + XU^ + Ux^) Q') 

+ K^oni{{x^u±u^x)QU^QU 



+ 1 



(2.9) 



where we have neglected inter alia terms proportional to 
(i.e. of second order in electromagnetic isospin breaking) 
that are disregarded throughout this work. The electro- 
magnetic LECs are renormalized in analogy to the strong 
LECs by 

K,=EiX + Kl{fi) , (2.10) 
with the coefficients Si given in Ref. [15]. 



contributions to the charged-pion and -kaon masses obey 
Dashen's theorem [31], 



(ZiM2)e^ = (ZiM|,)em + 0{e^mg) . 



(2.16) 



The Gell-Mann-Okubo (GMO) relation [32, 33] including 

isospin violation is given by 

2M|. + 2M|.o - 2M2 + M^o = 3M2 + 0{S^,p'') , (2.17) 

which is fulfilled in nature to a few percent accuracy. We 
define the GMO discrepancy according to 

2M| + 2MI-0 - 2M^ + M^o - 3M| 
^GMO = _ . (2.18) 



2.2 Meson masses and mixing at leading order 

By expanding the leading-order effective Lagrangian 
= C^l + one derives the LO meson masses. The 
flavor-neutral states 773 and % are mixed due to a differ- 
ence in the light-quark masses m„ — rrid ^ according 
to 



Bo ( TTs 



rriu + rud 



i=(m„ - rud) \ /tts 



V3 



2 XVsJ \-^{mu-md) l{mu + md + 'irns)j XVsJ ' 

(2.11) 

which can be diagonalized by the single rotation angle^ 



1 ^ /a/3 md - m„\ 

- arctan — — = —r 

2 \ 2 ms-m J 4 



md - ruu 



TTls — m 

(2.12) 

where for convenience and later use we defined the average 
light-quark mass rh and the isospin breaking parameter S, 



m ■ 



rUu + md 



S = md-mu 



(2.13) 



leading to the physical meson mass eigenstates 



We will frequently use the leading-order relations 

3M2 -M\ 
Bom,= \ '-° +0(.5^/), 



(2.19) 



as well as (2.15) to replace quark masses and Z by physical 
meson masses. 



2.3 Quark mass ratios 

The results of the previous section allow to form quark 
mass ratios determined entirely in terms of observable me- 
son masses [34], 



rrid 



Ml, -Ml + M^ 



0.55 , 



m« _ M|o +Ml-M^ 
md ~ M|o-M| + M2 



20.1 



(2.20) 



which are, however, subject to substantial higher-order 
corrections. The double ratio 



cos e sm e \ tt^ 
- sin e cos ej \ris 



(2.14) 



Mixing of the octet rjs and the singlet rjo to the observed 
mass eigenstates 77 and ry' is encoded in the strong NLO 
low-energy constant L7 [30] ; the heavy t]' is not a dynam- 
ical particle in this theory. 

For later use we define the charged-to-neutral meson 
mass differences 

AM^ = Ml - Ml, = 2F^Ze'' + 0{S^,p^) , (2.15) 
AMI = Ml- Mlo = 2F^Ze^ - B06 + Oip") , 

where M^, Mk denote the charged-pion and -kaon masses. 
Prom (2.15) it is obvious that the LO electromagnetic 



^ Note that this description of the ■k'^ti mixing via a single 
rotation angle is only valid at lowest chiral order [7]. 



ml -ml 



(2.21) 



nevertheless is particularly stable with respect to strong 
higher-order corrections [7], 



Ml Ml - Ml 



g^°^ (M»:-M^U <'"°K.^.O}. (2.22) 

Q is the major semi-axis of Leutwyler's ellipse [35] (ne- 
glecting a tiny term proportional to {rh/mgY), 



ma J \mdj 



(2.23) 



Furthermore, at leading order is also invariant under a 
shift in the quark masses of the form m„ — > m„ -|- a mdm.s 
{+ cyclic) [36]. 
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Fig. 3.1. LO diagrams for the charged and neutral decays. 



By using Dashcn's theorem in order to correct for the 
electromagnetic mass contributions, the doul)le ratio 
can be calculated at LO yielding a numerical value of 
Qd ~ 24.2, which corresponds to the LO quark mass ra- 
tios (2.20). However, Dashcn's theorem is subject to po- 
tentially large higher-order corrections and different mod- 
els (see e.g. Refs. [37-39]) yield a range 



1< 



(Mi 



(M2 - M^o) 



(2.24) 



which leads to a rather large uncertainty in the numerical 
value of Q of 20.6 < Q < 24.2 « Qd. 

Due to Sutherland's theorem [1], the dependence of 
the decay 77 — > Stt on the light-quark mass difference md — 
m„ is much less prone to be obscured by electromagnetic 
effects. As will be shown below, the decay amplitudes can 
be written in terms of Q^, and hence an accurate study of 
this decay can lead to an independent determination of the 
quark mass ratios. The re-evaluation of electromagnetic 
corrections in the present work allows for an increased 
precision in this determination. 



3 The decay amplitudes 

In this section we explain how to compute the decay ampli- 
tudes for both the charged and the neutral channel of 77 ^ 
37r decays at next-to-leading chiral order considering 
isospin breaking up to and including order e^{md — rUu), 

(ttV+tt- out|r?) = i{2wfS\P^,^-p„)Ac{s,t,u) , 

{3n° out|77) = i{2TT)^d\P^,, - Pr,)An{s, t, u) . (3.1) 

More technical details of the calculation can be found in 
Ref. [40]. 



3.1 Leading-order decay amplitudes 

At leading chiral order one only has to compute one tree 
graph each since the fields in the lowest-order Lagrangian 
are already diagonalized by use of the mixing angle e. 

The Feynman diagram describing the charged decay at 
lowest chiral order is shown in Fig. 3.1. The Mandelstam 



variables 



s = {Pv-P^of , t = {pn-Pn+f , u = iPv-PT^-f (3-2) 



are related by 



s + t + u = M^+M'^o + 2M^ = 3sg 



(3.3) 



where we defined the usual abbreviation sg. Inserting the 
mixing angle e, expanding to leading order in S, and sub- 
sequently replacing Bonis, Bom, and Z at leading chiral 
order by physical meson masses via (2.19) as well as Fq by 
the pion decay constant yields the well-known lowest- 
order result for the amplitude, 



A. 



LO 



1 



3V3F2 
(3s-4M2„)(3A/f2 



:i{s - + 2AMI 



- Ml, 



(3.4) 



An additional electromagnetic term of order €^{01^ — m„) 
cancels the pion mass difference implicitly included in Sq. 
The charged LO amplitude is completely proportional to 
md—rriu ^ l/Q^- It depends linearly on s, and by inserting 
Sq it explicitly displays the Adler zero at s = AM^a/i- 

The neutral decay is described at leading order by the 
diagram shown in Fig. 3.1, and in analogy to the charged 
decay we define the Mandelstam variables according to 

s= {Pn-P^of , t= {pr,-p^of , u= [p^-p^of , (3.5) 



which are related by 



s + t + u- 



Ml + ZMlo 



3s^ 



(3.6) 



The lowest-order amplitude contains neither derivatives 

nor electromagnetic terms and the consistent expansion 
in 5 and the replacements described above yield 



/iLO _ _ 



3(M^^-M^o)(3Af,^ + M^o) 



(3.7) 



The neutral LO amplitude also has an overall factor of 
"^d — but is just a constant. 



3.2 Contributions at next-to-leading order 

While for the charged process invariance under charge con- 
jugation implies that the amplitude Ac{s, t, u) is symmet- 
ric under the exchange of t and u. 



Ac{s,t,u) = Ac{s,u,t) , 



(3.8) 



the amplitude An{s,t,u) for the neutral decay has to 
be symmetric under exchange of all pions and thus all 
Mandelstam variables due to Bose symmetry. All calcula- 
tions of both ?7 — !■ Stt decay channels performed in ChPT 
so far [8, 13, 14] obey the following relation between the 
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charged and the neutral amplitude that utilizes isospin 
symmetry and the selection rule AI = 1, 

An{s, t, u) = Ac{s, t, u) + Ac{t, u, s) + Ac{u, s, t) . (3.9) 

For the LO amplitudes (3.4) and (3.7) this relation can 
be verified explicitly by use of (3.6). However, the rela- 
tion (3.9) is only valid at leading order in isospin breaking 
and cannot be used at 0{e^{md — m„)). This is most eas- 
ily seen by the fact that e.g. photon loop contributions do 
not respect this relation. Thus we have to calculate the 
neutral channel explicitly. 



3.2.1 General procedure 

In the following we will collect the contributing Fcynman 
diagrams up to next-to-leading chiral order p** for both 
decays. In order to obtain an unambiguous and explicit 
ordering of the isospin breaking parameters {(5, e^} with 
a minimal set of different meson masses we apply the fol- 
lowing procedure. 

1. We rewrite the physical meson masses induced by 
derivatives in terms of M^, M^o, ^ and AM"^ ~ Ze^ 
by using the GMO relation (2.17) and the charged-to- 
neutral mass differences (2.15). Higher chiral orders 
can be neglected everywhere except for the lowest- 
order diagrams, which we deal with explicitly in 
Sect. 3.2.2. 

2. We insert the mixing angle e (2.12) and expand every- 
thing in terms of 5. 

3. In order to disentangle the contributions to the various 
orders in isospin breaking ((5, e? , and (5e^), we expand 
the differences of certain charged and neutral kaon loop 
functions, both for the tadpoles Ak — A^o and for the 
two-point functions Jkk — Jk°k° (see (B.2), (B.4), 
(B.5) in Appendix B.l). We make use of the follow- 
ing expansion formulae, expressed in terms of 6 and 
AMI = Oie"^)- 



Ak - Ak^ = {AMI - BoS) 
BoSAM^ 



167r2 



167r2M|.o 



(3.10) 



T T W-AM^) 



^ifOKo(s)- — 



+ (.-4M^„)^ 

+o(e^<5^/). 



is) 



327r2M2.„ 



(3.11) 



We remark here in passing that kaon mass effects in 
loops as the above are the only source of terms of 
0{8'^) in — > 37r amplitudes. These are not of elec- 
tromagnetic origin and therefore disregarded. Their 
suppression with respect to the leading terms is never 
enhanced compared to the numerical size of e or 
AM'^/Mj.. 



Fig. 3.2. Diagrams contributing to yr^r/ mixing at NLO. The 
quotation marks serve as a reminder that the -k^ and the ri 
cannot be on-shell at the same time. 



4. We collect the terms with i5, e^, and Se^ and neglect 
all terms proportional to i5^, and of higher isospin 
breaking orders. 

5. We factorize l/Q^ ~ (5 where possible, using 



Bg5 = 



_ 3(M2 



Ml,){iM^^+Ml,) 



+ 0{p^) . (3.12) 



6. Finally, we replace left over quark masses and factors 
of Z by the LO meson masses M^, M^o, and AM^ 
via (2.19) and (2.15), and Fq by F„. 

This scheme yields a representation in terms of the isospin 
breaking parameters {l/Q^, e^, Z\M^} and the meson 
masses M^ and M^o . The latter choice is motivated by the 
fact the these masses appear as the asymptotic states in 
the processes under investigation, and feature naturally at 
next-to-leading order in tt^t] mixing, see Sect. 3.2.2. Note 
however that no mass substitutions take place inside loop 
functions, such that all cuts, even outside the physical re- 
gion, remain at their exact places. 



3.2.2 Tree diagrams and mixing 

The formalism to calculate matrix elements in the pres- 
ence of mixing, using the Lehmann-Symanzik-Zimmer- 
mann reduction formula, is explained in detail in Refs. [13, 
41]. Here we use an equivalent way to treat mixing, valid 
at one-loop order in ChPT, based on performing the cal- 
culation in a basis of states diagonalized at tree level. 

Working in terms of tree-level eigenstates, one has to 
account for tt^t? mixing effects beyond the leading order 
relation (2.12), (2.14) at chiral order p'^, induced by the 
mixing diagrams Fig. 3.2. Calculating the tt^t; self-energy 
matrix in the basis of tree-level eigenstates to one loop 
yields the off-diagonal element 



■^TT^riKP ) — >^TT"ri + ^iT°riP i ^7r"j) — 



, (3.13) 



written in a fashion that Z^o^ can be interpreted as the 
off-diagonal analogy to the wave-function renormalization 
factors. E.g. in the charged decay channel, NLO mixing 
in the outgoing 7r° leg can be written in terms of the LO 
amplitude ^J^^-r/Tr+Tr- '^ith subsequent mixing of the out- 
going 7r° due to the diagrams of Fig. 3.2, 



A: 



LO SnOniM^o) ^ ^LO 



■ £4 



1 [Ko+Ml 

~ M\-MS^M 2 



(3.14) 
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a) Tree diagram 



b) Counterterm diagram 



Fig. 3.3. NLO trcc-typo diagrams. The double solid line de- 
notes the full propagator, cf. Fig. 3.4. 



(see e.g. Ref. [42]). In analogy the mixing of the ingoing ry 
with a 7r° is given by 



M2 . 
-'"77 



A: 



LO 



£4 



k 2 

(3.15) 

Since both Z^Oj^ and £4 are of chiral order and of first 
order in isospin breaking, all quark mass replacements in 
the corresponding amplitudes can be done at order ne- 
glecting terms proportional to 6e'^. The expressions for the 
corresponding LO amplitudes as well as Z^o^ and 64 are 
displayed in Appendix A.l. Taking into account the mix- 
ing of the ingoing 77 and all outgoing tt'^ , the NLO mixing 
contributions for the charged and neutral channel can be 
summarized in the form 



^mix 



3s 



M2 



Ml. „ 



£4 



£4 



(3.16) 



With the mixing phenomenon dealt with, we can 
continue by considering the tree diagram displayed in 
Fig. 3.3a) consisting of a vertex of order p^ and full prop- 
agators up to order p^. The NLO propagator that is de- 
picted in Fig. 3.4 contains the LO propagator and proper 
self-energy insertions —iSip^) of order p** from meson tad- 
poles, from a graph with a vertex from C'^^\ and from 
photon loops for the charged pions. The sclf-cncrgics at 
chiral order p'^ allow to compute both the masses at NLO 
needed for renormalization corrections to the substitutions 
of quark masses by meson masses, and the wave-function 
renormalization factors 



Za 



dSaip^) 



= l + AZ, + 0{p^). 

dp'' p^=Mi J 

(3.17) 

Specifically, S-^a and refer to the diagonal elements of 
the self-energy matrix in tree- level eigenstates. By defining 
the following renormalization corrections 



M2 - m2 = 



4Bo 

o I?2 /7„2 



{nris -rh)+ Am , 
2F^Ze^ + Az , 



3M2 + M^o rris 



m 



{1 + ^0} 



2™ 

F^=Fo{1 + AfJ , 



(3.18) 



the corresponding NLO renormalization corrections to the 
LO amplitudes read 



3M,f + Mio 



-Am~2Az] , 



^[{2.s-^MI,){2Ap^-Aq) 



^ren _ 



^/3(3M2 -HM^o) r „ 
-Am] ■ 



(3.19) 



The formulae for the Z factors and the renormalization 
corrections defined above are shown explicitly in Ap- 
pendix A.l. Since the value for the physical pion decay 
constant F^^ = 92.2 Me V [43] was extracted from decays 
of charged pions with electromagnetic corrections already 
taken care of (sec also Rcf. [44] ) , we use the charged-pion 
decay constant in the absence of electromagnetism. The 
relation to the bare Goldstone boson decay constant Fq is 
hence as given in Ref. [7]. 

Besides the renormalization effects discussed above, 
both strong and electromagnetic NLO low-energy con- 
stants enter the calculation via the diagram Fig. 3.3b) 

with a vertex from = C^^^ + £om- We denote the re- 
sulting contributions by and , respectivcily. For the 
charged channel one finds a contribution proportional to 
is quadratic in the Mandelstam variables, terms linear 
in s, and constant terms. However, the contribution for 
the neutral channel does not depend on L3, L4, and Ki 
to (i.e. terms with traces containing only derivatives). 
Therefore it is a constant like the LO amplitude as well 
as other tree contributions. 

Combining all parts discussed in this section, the tree 
diagram contributions up to next-to-leading chiral order 
then read 



1 



AZ„ 



+ 



1 + 



AZ, 



+ ■ 



AZ^ 
2 

■^AZ^o 



AZ^ \ A. 



iLO 



jrcn I ^mix , Act 



iLO 



2 2 

, ^ren , ^mix , ^ct 



(3.20) 



In ^'"/"j, any explicit dependence on Z/4 and Lq cancels. 
The remaining strong LECs are replaced in terms of ob- 
servables in the spirit of Ref. [8]. Besides the correction 
^GMO to the Gell-Mann-Okubo relation (2.18), we define 
the correction Ap to the ratio of the kaon and the pion 
decay constants at NLO without isospin breaking as given 
in Ref. [7] according to 



Fk 



l + Ap. 



(3.21) 



The explicit forms of these two quantities are given in 
Appendix A.l. Since both A^^J'^ depend only on the com- 
bination 2L7 -|- Lg; both L7 and Lg can be eliminated si- 
multaneously in terms of Aqmo , while L5 can be replaced 
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+ + . + 

Fig. 3.4. Full propagator up to 0{p^). Wiggly lines denote photons. 




a) Tadpole diagrams b) Rcscattcring diagrams 
Fig. 3.5. Meson loop contributions. 

by Ap. The explicit formulae for A^J^'^ and A^^'^^ follow- 
ing from these replacements are given in Appendices A. 2, 
A. 3. In this form, the neutral amplitude has no explicit 
dependence on any strong LEG, while the charged one 
only contains a term proportional to L3. Furthermore, as 
a by-product, in both any explicit dependence on the elec- 
tromagnetic constants and Ks cancels. 

3.2.3 Meson loops 

At chiral order p'^ hadronic loops can occur in three places. 
In addition to the tadpole contributions to the NLO prop- 
agator that have been discussed already, there are also 
diagrams with a meson tadpole directly at the LO ver- 
tex and diagrams with rescattering intermediate mesons 
as depicted in Fig. 3.5. 

The results for the tadpole contributions A^^ and 
A^n'^ arc given in Appendices A. 2. A. 3. Again the charged 
contribution contains terms linear in s and constant terms, 
whereas the neutral contribution is constant. 

The rescattering of intermediate mesons, in particu- 
lar / = TTTT final-state interaction, is very important and 
gives rise to about half of the total NLO corrections for the 
decay 77 ^ Stt [8]. As pointed out before the charged decay 
is symmetric under t <-> u and the neutral decay is sym- 
metric binder exchange of all Mandclstam variables. While 
for the neutral decay it turns out that the contributions in 
each channel depend on the corresponding variable only, 
in case of the charged decay both the t and the u chan- 
nel contributions depend on at least two variables. Since 
there are no virtual-photon loops for the neutral decay, 
these are the only dependences on Mandelstam variables. 

Defining the functions A'l, A'^" and A^^ as explained 
in Appendices A. 2, A. 3, the rescattering contributions for 
both decays read 

_^scatt ^ At!'"is) + ^^'"(t) + ^^'"(m) . (3.22) 

3.2.4 Photon loops 

For the charged decay at chiral order p^, photon loops oc- 
cur. These typically entail infrared (IR) divergences, which 




a) Vertex correction diagrams b) Triangle diagram 
Fig. 3.6. Virtual-photon loop contributions. 

are cancelled in the standard manner by the inclusion 
of real-photon radiation (bremsstrahlung) on the cross- 
section level. We keep track of the infrared divergences 
by introducing a small photon mass niy-, alternatively one 
might use dimensional regularization also in the infrared. 

In addition to the self-energy contributions to the NLO 
propagators of the charged pions. there are also vertex cor- 
rections as depicted in Fig. 3.6a), and the triangle diagram 
displayed in Fig. 3.6b), describing the exchange of a vir- 
tual photon in the final state between the charged pions. 

The photon loop contribution U^{p'^) to the self-energy 
of the charged pions contributes both to the (IR-finite) 
mass and the (IR-divergent) wave-function renormaliza- 
tion. The two vertex correction diagrams turn out to yield 
the same IR-finite contribution A'^^ and the triangle dia- 
gram leads to the contribution A^'^^ . Both are given ex- 
plicitly in Appendix A. 2. The latter one contains the tri- 
angle loop fimction G{s) (B.7) which has some interesting 
features: both the real and the imaginary part are IR- 
divergent; while the infrared divergence in the real part is 
cancelled against bremsstrahlung contributions, the imag- 
inary part can be resummed in the (divergent) Coulomb 
phase. Furthermore G{s) contains a kinematical singular- 
ity at threshold s = 4M^, the Coulomb pole. Its contri- 
bution to the complete amplitude is given by 

^pole^^LO ^g2l+^ (3.23) 
c c 16C7 ^ ' 

where tr = a/1 — 4M^/s. Note that the prefactor of the 
triangle loop function, i.e. the Coulomb pole, is propor- 
tional to the LO amplitude (3.4). 

3.2.5 Real-photon radiation 

Infrared divergences due to virtual-photon corrections are 
canceled by including real-photon radiation; see Fig. 3.7. 
We only include real-photon radiation in the soft-photon 
approximation that amounts to neglecting the photon mo- 
mentum in the overall energy and momentum conserva- 
tion. The results for general n-body decays with one addi- 
tional real photon of maximum energy i?max radiated can 
be found, e.g., in Ref. [45] (also using a finite photon mass 
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Fig. 3.7. Real-photon radiation diagram. 



as a regulator); radiative corrections for the largely anal- 
ogous decay Kl tt+tt^tt^ in the framework of ChPT 
are discussed in Ref . [46] . Adapting the results of Ref. [45] 
to T] ^ 7r+7r~7r°7 yields an effective contribution to the 
amplitude squared of the form 



\Ae 



47r2 



In 



4£;2. 



2cr 



In 



1 



F{s,t, 



The maximum photon energy E„ 
is given by 



E„ 



minjE'kin, -Bcut} > -E^kh 



1-a 

(3.24) 

(in the r] rest frame) 



- (M^o + Vi)' 



(3.25) 

where -Ekin is the maximal kinematical limit, and the value 
of the photon cutoff energy i?cut is set to a typical detector 
resolution. The kinematical constraint leads to a logarith- 
mic divergence at the upper limit s = (M^ — M^o)^, sig- 
nalling that very close to the boundary of phase space the 
0{a) approximation becomes unreliable. F{s,t,u) repre- 
sents contributions that are finite in the limit nij 0. 
They are given explicitly as 

F{s,t,u 



uj{z) 




zv -\- uj(z) 



zv 1 
■ m 



2a>(z) \ + zv — {jj(z) 
y^P^ + 2ziy + z'^{4p^ + 



Ml 



Ml 



Mr,^s-mi 

M^ + s- Ml, 



(3.26) 



where A(a, 6, c) — + + — 2{ab + bc+ ca) is the Kallen 
function. 

Computing the amplitude for 7r+7r~7r° squared 

and collecting the IR-divergent parts consistently with re- 
spect to chiral order and power counting in we obtain 



' 47r2 Ml 



1 In 

2a 1 - C7 



(3.27) 



Comparing (3.24) and (3.27) shows that the inclusion of 

bremsstrahlung cancels the IR-divergences. 

We have been intentionally vague above about the am- 
plitude Ac to be inserted in (3.24), (3.27). As a matter of 
principle, our calculation is only fully consistent for the 



lowest-order amplitude A}^'^ , to be multiplied with electro- 
magnetic corrections factors. However, as is conventionally 
done in ChPT calculations, we only chirally expand the 
amplitudes to a specific chiral order, and we do not re- 
expand their squares. This way, even when disregarding 
terms of order (which are tiny), we include interference 
terms of radiative corrections with strong loop corrections. 
As a consequence, in order to achieve cancellation of in- 
frared divergences, the bremsstrahlung terms in (3.24) are 
included with a prefactor 



^^°(^^° + Re{^^^°|e=o}) 



(3.28) 



3.2.6 Subtraction of universal corrections 



The kinematical singularities in the radiative corrections 
both at s = 4M2 (3.23) and at s = (M^ - M^af (3.24), 
(3.25) are part of the universal soft-photon corrections 
that can even be rcsummed to all orders in the fine struc- 
ture constant [45, 47, 48] . In order to perform a meaningful 
fit of the Dahtz plot distribution (see Sect. 4.3), these uni- 
versal soft-photon corrections are usually already applied 
in the analysis of the experimental data. We account for 
this by subtracting the corresponding parts of our am- 
plitude (squared) before predicting experimental observ- 
ables, in the following way. 

1. The both infrared and kinematically divergent imag- 
inary part of the photon triangle loop function G{s) 
(B.7) can be rcsummed to the (divergent) Coulomb 
phase. As an overall phase factor is unobservable, 
Im{G(s)} is omitted. 

2. We subtract the Coulomb pole, whose leading approx- 
imation in our calculation is given in (3.23). 

3. We subtract the Ini^max singularity of (3.24) in the 
form 



1 2a 



(3.29) 



which is the 0{e^) approximation of the resummed 
correction factor, see Refs. [45,48]. 

In this manner, we obtain a squared amplitude free of 
kinematical singularities. 

For the purpose of illustration of the size of the various 
corrections, we find it useful to plot (real and imaginary 
parts of) amplitudes in Sect. 4.2, which, as illustrated in 
the previous section, can in principle not be made infrared 
finite in a consistent way. We remedy this problem by 
hand, using the replacement 



(3.30) 



in the decay amplitude, compare (3.24), (3.27), (3.29). 
This does not take into account the finite contributions 
F{s, t, u), which can only be added on the level of squared 
matrix elements. Furthermore, in Sect. 4.2 we retain the 
threshold divergences (Coulomb pole and phase) for illus- 
trative reasons. 
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3.3 Total decay amplitudes 

In terms of the contributions defined before and by ap- 
plying the replacement rules given above the decay ampli- 
tudes finally can be written as 



SP 

where = denotes the evaluation at the soft-pion point 
and the Kl are to be evaluated at the scale /i = Mk- 
The additional terms of 0{e'^6) (denoted by DKM), for 
comparison, are found to be 



^tree 
^trec 



A tad 



^scatt _j_ -|- ^ 

■ . (3.31) 



Note that in case of the neutral decay the whole kinemati- 
cal dependence is contained in the scattering contribution, 
which takes the symmetric form (3.22). 

We have checked both the charged and the neu- 
tral amplitude in several ways. They are finite and 
renormalization-scale independent. In addition, we have 
explicitly verified that both amplitudes reduce to the re- 
sults of Ref. [7] at 0{5) and of Ref. [14] at 0{e^). 

We have pointed out before that the total amplitudes 
(3.31) do not obey the AI = 1 selection rule (3.9). The 
terms of 0{6e^) violating this rule are, however, difficult to 
write down or even define consistently, as the pion mass 
difference affects even the relation between the Mandel- 
stam variables in the two channels (due to sg ^ Sq). Just 
for the purpose of illustration, we show here the part of 
the deviation proportional to the various low-energy con- 
stants, using s + t + u = 3so : 



A 



LEG 



A^^""^ {s, t, u) ~ A^^^ {t, u, s) - yl^^^ (u, s, t) 



3M^ + M^o 

1{2K; - Kl){Ml + Mlo) - K^,M^ 



(3.32) 



As a final analytic result, we want to comment on the 
soft-pion theorem [1] for the electromagnetic corrections. 
As shown in Ref. [14], the corrections of 0{e^mq) (where 

m„ = rud is assumed) at the kinematical point s = 3sq^", 
t = u = (the soft-pion point) scale as e^m, not as e^nis- 
Explicitly, we find for the electromagnetic amplitudes in 
this approximation (denoted by BKW) 



A. 



BKW SP 



-F^{ 1 - \/2arctan —= 

4e2 



+ 



-i2Kl - Kl) 



Kl 



Kl + Kl + K{o 

(3.33) 



A. 



BKWSP ^^"/ ^^|^(^i_^,,etani=) 



+ 



V3F2 1 16^2^2 
4e2 



Kl 



Kl 



Kl + K\, 
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(3.34) 



A, 



dkmSp 2Bo5\AMl 
^^^F^\ M2 



A, 



GMO 



At 



F2 



+ 



m: 



167r2F2 V3 Mk 2 3 

7 1 „ . 

H ;= arctan — ;= + 2ni 

4V2 V2 



2 r 



e 
8^ 



(l + 21n^+.^)ln^-ln^^ 



TT^ 9 TT . , 

1 V - i + — - 

3 4 4 Mk 



V 

4 3 



M^ 



4e2 
~3~ 



7 



K{ + Kl, - -{2KI - Kl) + Kl + Kl 



-KE + 2Klo + 2,K 



A. 



DKM SP 2Bo 6 f AM? 



+ e 



-),(3.35) 



V^F^ I M2 



M2 



3^^^ 167r2f2 



In- 



M^ 



K 



^,4^1 1 

-|- 2 In — 1-3-1 ;= arctan — ;= 

3 4\/2 V2 



4e^ 
~3~ 



^+K[ + K^, --{2Kl-Kl)+ Kl 



Kl-Kl + 5Klo + 6Kl, 



O 



(3.36) 



again at = Mk- The comparison of (3.33)-(3.36) 
demonstrates explicitly that terms of 0{e^6) are relatively 
suppressed only by (m^ — m,y)/{'md + ni-a) ~ 1/3 and 
not by another small isospin- violating parameter. Further- 
more, only the terms in (3.35), (3.36) show chiral loga- 
rithms oc InM^. The chiral logarithms squared in (3.35) 
are due to the photon triangle loop function G{s), ex- 
panded around the soft-pion point. 



4 Numerical results 

The numerical input that was used to obtain the following 
results is collected in Appendix C. 



4.1 Kinematical bounds and cusps 

In Fig. 4.1 we plot the boundaries of the Dalitz plot for 
both decay channels. In addition to the charged and the 
neutral decay phase spaces we have also drawn the allowed 
phase space for the decay 77 Tr+Tr^Tr*^ using an average 
pion mass A?2 — {M'^o + 2M2)/3. This approximation is 
used frequently in the literature, since it is, strictly speak- 
ing, not possible to consistently account for different pion 
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s [GeV^] 

Fig. 4.1. Kinematical bounds of the Dalitz plots for the decays 
7] —>■ TT+TT^Tr" (full/red), rj — > Stt" (dashed/blue), and rj 
7r^7r~7r° with an average pion mass (dotted/magenta, see 
main text for details). Furthermore shown are the lines with 
t = u for both the charged and the neutral decay as well as 
the 4:M^ thresholds in all three kinematical channels (dashed- 
dotted/black). 

masses (which feature in 77 Stt at 0((5e^)) at leading 
order in isospin breaking. Deviations at the border of the 
Dalitz plot are clearly visible. 

Recently, there has been a strong renewed interest in 
threshold cusp phenomena [49] (for an overview of var- 
ious examples, see also Ref. [50]) since Cabibbo [19,20] 
pointed out that the large-statistics experimental data 
on the decay K+ 7r°7r°7r+ collected by the NA48/2 
Collaboration at CERN may allow for a very precise de- 
termination of the pion-pion scattering length combina- 
tion oo — 02 [21]. These scattering lengths multiply the 
strength of the cusp in the tt^tt" invariant mass distri- 
bution at the tt'^tt^ threshold; they are, at the same 
time, predicted theoretically with tremendous precision, 
Uq — a2 = 0.265 ± 0.004 [51, 52], and represent a core test 
of ChPT and our picture of chiral symmetry breaking. The 
tailor-made theoretical framework to perform a precision 
analysis of these cusps is non-relativistic effective field the- 
ory [23-25] that has been exploited up to two loops and 
including radiative corrections. 

A similar cusp phenomenon also occurs in the decay 
77 Stt*', although, as in the analogous decay 
377° [53], it is less pronounced than in A'+ — > 7r"7r"7r+. 
From an experimental point of view 77 —> Stt^ decays have 
therefore not been competitive for a pion-pion scattering 
length determination. But as new experiments like WASA- 
at-COSY [54, 55] or Crystal Ball at MAMI-B [56-59] will 
provide more high-statistics data in the near future, an 
investigation of the cusps also in 77 — s- Stt*^ might become 



feasible. It has already been explored theoretically (along 
with Kl ~» 37r") in the non-relativistic framework [24,25] 
(for a qualitative study using a relativistic field theory, 
see Ref. [60], for an assessment in the framework of uni- 
tarized ChPT, see Ref. [61]). As the present calculation 
is performed in chiral perturbation theory, and hence in- 
cludes the rescattering effect leading to the cusp only at 
leading order in the quark mass expansion, it is not suited 
to serve for a precision extraction of pion-pion scattering 
lengths, but will rather illustrate the phenomenon. It is in 
a sense dual to the non-relativistic calculation as it aims 
at predicting electromagnetic effects in those parts of the 
amplitude that are merely parameterized in the latter. 

In Fig. 4.1 we show the cusp lines for the energy thresh- 
old of the production of two charged pions in all three kine- 
matical channels s, t, and w as a vertical, horizontal, and 
diagonal dashed line. Furthermore, for illustration pur- 
poses we plot the amplitudes along the lines t = u in the 
following section; these lines are also drawn in Fig. 4.1. 

4.2 Decay amplitudes 

In this paper we focus on the electromagnetic contribu- 
tions to 77 ^ 37r decays and do not aim at a particularly 
reliable representation of the purely strong amplitude. It 
is well known that one has to go beyond one-loop order 
to obtain a valid representation of the latter [9-13], while 
the 0{S) part in our calculation corresponds precisely to 
the one-loop representation of Ref. [8]. This serves as a 
useful reference point to quantify the electromagnetic cor- 
rections. In the following, we only consider uncertainties 
in precisely these electromagnetic contributions and dis- 
regard higher-order hadronic corrections. Since the elec- 
tromagnetic low-energy constants are not very well 
known, we regard their input (and not unknown correc- 
tions of higher order in the chiral expansion) as the domi- 
nant source of uncertainty; see Appendix C for a descrip- 
tion of how we vary the K^. For some numbers we are 
forced to deviate from our standard procedure to estimate 
the errors, those are marked by an asterisk in the follow- 
ing; the details are also discussed in Appendix C. All error 
bands in this section refer to this variation. 

In the following, we compare the results for the one- 
loop amplitudes of 0{6) (henceforth denoted by GL) [8], 
those with effects of 0(6^) added (BKW) [14], and the 
complete results of the present investigation, up-to-and- 
including effects of 0{de'^) (DKM). The GL and BKW 
amplitudes are also evaluated with our prescription of the 
isospin limit for the pion mass (as well as with our choice 
of numerical input. Appendix C) in order to facilitate com- 
parison with the higher-order corrections. 

In Fig. 4.2 we separately display the real and the imag- 
inary parts of the charged amplitudes GL, BKW, DKM 
along the line t = u. The infrared divergences in the am- 
plitude are cured by hand according to (3.30); the kine- 
matical singularities at s = 4M^ are retained here for 
illustration. By our choice of the isospin limit for the pion 
mass, the threshold cusp in the GL amplitude is artificially 
removed from the physical threshold energy to s = 4M^o ■ 
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Fig. 4.2. Real and imaginary parts of the charged ampli- 
tudes GL (dashed/black), BKW (dot-dashed/blue), and DKM 
(full/red) for t = u. The inserts show the region close to 
the two-pion thresholds. The line widths in the real part in- 
dicate the error bands due to a variation of electromagnetic 
low-energy constants. The vertical lines show the limits of the 
physical region. 



According to (3.4) the leading-order charged decay am- 
plitude is linear in s, which still dominates the energy 
dependence at NLO. The BKW contributions are purely 
real for both decay channels, as no pion rescattering di- 
agrams contribute at that particular order and thus the 
imaginary parts of GL and BKW coincide with each other 
for both decays. Furthermore, all imaginary parts are in- 
dependent of any low-energy constants and are therefore 
plotted without an error range. 

Since it is hard to identify cusps or the Coulomb pole 
in plots over the full kinematically allowed range. Fig. 4.2 
also contains inserts where we show the amplitudes close 
to the two-pion thresholds. Here the expected features in 
the amplitudes are clearly visible: the tt^tt^ cusp at 4M^o 




0.06 0.08 0.10 0.12 0.14 0.16 0.18 
s [GeV^] 




\ 



0.06 0.08 0.10 0.12 0.14 0.16 0.18 
s [GeV^] 

Fig. 4.3. Real and imaginary parts of the neutral ampli- 
tudes GL (dashed/black), BKW (dot-dashed/blue), and DKM 
(full/red) for t — u. The hatched regions denote the error bands 
due to a variation of electromagnetic low-energy constants. The 
vertical lines show the limits of the physical region. 



outside the physical region, as well as the Coulomb pole 
and phase divergence at the tt^tt" threshold. 

The numerical results for the neutral decay amplitudes 
GL, BKW, and DKM are shown in an analogous way in 
Fig. 4.3. The vertical dotted lines indicate the physical re- 
gion in s for t — u for the neutral decay, hence the lower 
bound is now at the tt^tt" threshold. The leading decay 
amplitude for — > 37r*' is constant, see (3.7), and also at 
NLO the dependence on s is weak. Figure 4.3 shows that 
the size of the new contributions in DKM is comparable to 
or even larger than those in BKW. We find the expected 
cusp at the energy of the tt+tt^ threshold inside the phys- 
ical region, although the overall variation is very small. As 
for the comparable decay Kl 37r°, the strength of the 
cusp might be rather small for a very precise determina- 
tion and an extraction of tttt scattering lengths thereby. 
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We wish to point out, however, that the cusp strength 
is underestimated here. As we mentioned before it is pro- 
portional to the combination of scattering lengths oq — 02, 
and the present ChPT calculation corresponds to the 
LO values (oq - 02)^° = 9M2/(327rF2) = 0.204 [62]^ 
while the matching of the Roy equations solution to 
the next-to-ncxt-to-lcading chiral representation yields 



(ao-a2) 



NNLO 



0.265 ±0.004 [51,52], increasing the cusp 



strength considerably. Furthermore, it has been noted be- 
fore that due to the relative strength of the decays into 
charged and neutral final states, the cusp is much less 
pronounced in 77 ^ Stt^ than in tt'^tt'^'it'^ , much 

as in Kl —> Stt" [24]. For these two reasons, we refrain 
from displaying the decay spectra dF/ds, which would be 
the preferred observable for an extraction of tttt scattering 
lengths. 

From Figs. 4.2 and 4.3 we can conclude that the rela- 
tive sizes of the electromagnetic corrections in the BKW 
and DKM amplitudes with respect to the strong result 
GL are of comparable size, and the estimate of the elec- 
tromagnetic effects, neglecting m„ ^ nid therein, is not 
a very accurate representation. However, the conclusion 
that the overall electromagnetic contributions remain very 
small is still valid. 



4.3 Dalitz plot parameters 

The Dalitz plot distribution for the charged decay channel 
is described in terms of the symmetrized coordinates 



2M^g, 



3To 3[{M, - M^or - s] 

where = Tq + T+ + T_ = - 2M^ - M^o and the 
Tj are the kinetic energies of the respective pions in the 77 

rest frame. For the neutral decay channel it is convenient 
to use one fully symmetrized coordinate 



i=l 



Qn 



x^ + y\ (4.2) 



with Qr, 



Mr, 



3Mj^o , in order to reflect the symmetry 



in all Mandelstam variables. 

The squared absolute values of the two decay ampli- 
tudes are expanded around the center of the corresponding 
Dalitz plot according to the standard parameterizations 

\Ac{x,y)\'' = |AAe|2{l + ay + 5j/2 + dx2 

+ fy^ + gx^y + ...} , 
\An{x,y)\'' = \Afnmi + 2az + ...} , (4.3) 



^ To be more precise, tiic cusp strength in our calculation is 
determined by the leading-order scattering length for tt^tt" 
tt^tt'' including isospin breaking, which is given by (00 — 02)'^'^ x 
{1 -I- (M^ — M^o)/(3M^)}, increasing the isospin-symmetric 
value by about 2%. 



in order to obtain the Dalitz slope parameters. For the 
charged channel odd terms in x are forbidden due to 
charge conjugation symmetry. The cubic parameters have 
only been measured by the KLOE Collaboration [63] with 
only / found to be different from zero with statistical sig- 
nificance. (For older experiments, see Refs. [64-66].) 

Note that at second order in isospin breaking one has 
to take care of the precise definition of the center of the 
Dalitz plot. For the charged decay the point s = t = u = 
Sq does not completely coincide with the point x = y = 
(where all kinetic energies are equal) due to the chargcd- 
to-neutral pion mass difference, see Ref. [67]. Indeed, s ~ 



t 



c/n 



corresponds to a; = y = for the neutral 
decay, but x = and y = 0.0514 for the charged decay. 

We fit the distribution functions (4.3) to discretized 
grids (of roughly 200 x 200 points) of squared amplitudes 
over the whole physical region. We use uniform weighting 
of the "data" points. In order to quantify the fit quality for 
the amplitudes GL, BKW, DKM, we normalize a fictitious 
to 1 for the GL amplitudes and only regard the rela- 
tive changes, with the electromagnetic effects switched on 
successively. As explained above, in order to estimate the 
errors in the corrections to the various Dalitz plot param- 
eters, we vary the electromagnetic low-energy constants 
Ki as described in Appendix C. The resulting errors are 
always larger than the fit errors on the extracted param- 
eters, hence we can neglect the latter. 

The results for the normalization and the Dalitz plot 
parameters a—g of the charged decay channel arc shown in 
Table 4.1. We display the fit results for the GL, BKW, and 
DKM amplitudes. Note again that the universal radiative 
corrections producing kinematical singularities have been 
subtracted in the DKM amplitude according to the pre- 
scription given in Sect. 3.2.6. 

Table 4.1 shows that, in general, electromagnetic cor- 
rections affect the Dalitz plot parameters at the percent 
level. The normalization tends to get reduced compared to 
the purely strong amplitude, while the various slope pa- 
rameters a — g are slightly increased. The relative shifts in 
d and g arc more sizeable than in particidar in a and b for 
the reason that the strong contribution to the a;-dependent 
Dalitz plot parameters is suppressed to next-to-leading or- 
der in the chiral expansion (as is obvious from the purely 
s-dependent tree-level amplitude (3.4)). While the over- 
all effects are still very small, we note that, throughout, 
the corrections of 0{e^) do not represent a valid estimate 
of the dominant electromagnetic corrections, as those of 
0{5e^) are of the same order of magnitude - sometimes 
with the same sign; sometimes both effects tend to cancel. 
W; quote the corrections both as shifts of absolute size 
and as relative shifts in percent, compared to the strong 
one-loop result [8]. Regarding the fact that the strong 
corrections at 0{p^) are potentially substantial [13], the 
comparison of both may be taken as another indication 
of higher-order terms in the electromagnetic corrections. 
The quality of the polynomial Dalitz plot fit is comparable 
throughout. 

As explained before, the uncertainties quoted in Ta- 
ble 4.1 are obtained by varying the electromagnetic low- 
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Table 4.1. Normalization and Dalitz plot parameters for rj 7r"'"7r~7r'' for the GL [8], BKW [14], and DKM (present work) 
amplitudes. All electromagnetic corrections are given as shifts (both absolute and relative, in percent) with respect to the strong 
result (GL). For details, see main text. For a discussion of the errors, see Appendix C. 







a 


b 




GL 


0.0325 


-1.279 


0.396 




BKW 


-0.0004 ± 0.0003 


-0.008 ± 0.001 


+0.006 ± 0.001 






= (-1.1 ± 0.9)% 


= (+0.6 ± 0.1)% 


= (+1.4 ± 0.2)% 




DKM 


-0.0008 ± 0.0002* 


-0.009 ± 0.005 


+0.006 ± 0.003 






= (-2.4 ± 0.7*)% 


= (+0.7 ± 0.4)% 


= (+1.5 ± 0.7)% 






d 


/ 


g 


xVndf 


GL 


0.0744 


0.0126 


-0.0586 


= 1 


BKW 


+0.0011 ± 0.0004 


-0.0003 ± 0.0001 


-0.0010 ± 0.0003 


1.03 




= (+1.5 ± 0.5)% 


= (-2.2 ± 0.4)% 


= (+1.7 ± 0.6)% 




DKM 


+0.0033 ± 0.0003* 


+0.0001 ± 0.0001 


-0.0038 ± 0.0009* 


1.63 




= (+4.4 ± 0.4*)% 


= (+0.5 ± 0.6)% 


= (+6.4 ± 1.5*)% 





energy constants Ki. In most cases, the errors given are 
of the same order as the central shifts. 

A few comments concerning the cubic Dalitz plot pa- 
rameters / and g are in order. In general, / is highly corre- 
lated with a in the fits. In the electromagnetic corrections, 
wc observe that the influence of the low-energy constants 
on / and g is rather small, which is obvious from the 
fact that the counterterms only lead to terms constant 
and linear in s. We therefore do not consider the electro- 
magnetic shifts thus obtained very reliable. In particular, 
all errors due to the in the parameters d, f, g are 
only indirectly induced by the variation in the normaliza- 
tion and by the implicit inclusion of higher-order effects 
in the squared amplitude. Note that, without subtracting 
the kinematical singularities due to universal corrections 
(see Sect. 3.2.6), in particular the fit results for / become 
nonsensical. 

We now turn to the neutral channel ry Stt^. The 
results for the normalization and the Dalitz plot parame- 
ter a are collected in Table 4.2. Here, in addition to the 
results using the GL, BKW, and DKM amplitudes, we 

also discuss a variant of our result that takes into account 
that the DKM amplitude displays features incompatible 
with a simple polynomial fit: the cusps at the charged- 
pion thresholds. By DKM(cusp) we denote a fit to the 
part of the Dalitz plot with z < zq, zq = 0.598 chosen 
such that the border region from the cusp outward is ex- 
cluded.^ As for the charged decay channel, the normal- 
ization is reduced by electromagnetic corrections by a few 
percent. Corrections of 0{5e'^) are even bigger than those 
of 0{e'^). We note that the cusp effect leads to the single 
biggest modification of any Dalitz plot parameter: trying 
to fit the cusp with the polynomial distribution function 
reduces a by 4% (compare DKM to BKW in Table 4.2), 

^ Note that zo corresponds to the minimum value of z where 
the cusps occur, as these show up for constant s/t/u, not for 
constant z; see also the discussion in Ref. [61]. 



Table 4.2. Normalization and Dalitz plot parameters for 
ri 37r° for the GL [8], BKW [14], and DKM (present work) 
amplitudes. Results for the latter arc also shown with the fit 
of the Dalitz plot region restricted as to exclude the cusps at 
the charged-pion thresholds (DKM(cusp)). All electromagnetic 
corrections are given as shifts (both absolute and relative, in 
percent) with respect to the strong result (GL). For details, 
see main text. 





iA/;i' 


10^ X a 


xVndf 


GL 


0.269 


1.27 


= 1 


BKW 


-0.003 ± 0.002 


+0.05 ± 0.01 


0.99 




= (-1.1 ± 0.9)% 


= (+3.7 ± 0.5)% 




DKM 


-0.009 ± 0.005 


-0.002 ± 0.01 


6.20 




= (-3.3 ± 1.8)% 


= (-0.2 ± 1.0)% 




DKM (cusp) 


-0.009 ± 0.005 


+0.06 ± 0.01 


0.35 




= (-3.3 ± 1.8)% 


= (+5.0 ± 1.1)% 





while excluding the cusp region increases it again by more 
than 5%. The latter shift is in qualitative agreement with 
the finding in Ref. [16], where a is determined just from 
the curvature at the center of the Dalitz plot. The signifl- 
cance of this non-analytic structure is also reflected in the 
fit quality as quantified by the x^/ ndf values quoted in 
Table 4.2: with the cusp included, the fit becomes worse 
by a factor of 6 (DKM), while excluding it makes it even 
better than the fit of the GL distribution (DKM (cusp)). 
Taking into account that the cusp strength is underesti- 
mated by about 30%, see the discussion in Sect. 4.2, these 
numbers should be scaled accordingly. On the other hand, 
the cusp effect is by far too weak to contribute significantly 
to an explanation of the long-standing sign discrepancy for 
a between ChPT calculations [8, 13] and experimental de- 
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1.03 



1.025 



Table 4.4. Branching ratios r = FnlTc- The corrections for 
the BKW and DKM amplitudes are expressed relative to the 
strong GL result. DKM(uc) denotes the variant of the ampli- 
tude without the universal photon corrections subtracted. 



1.02 



1.015 



1.01 







1.442 




( 


-0.1 ± 1.2)% 




( 


-1.4 ± 1.8)% 


^^DKM(uc) 


( 


-2.3 ± 1.8)% 



1.005 



-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 
X 

Fig. 4.4. r\ Stt" Dalitz plot distribution corresponding to 
the full one- loop amplitude (DKM), including the cusp struc- 
tures at the 7r"'"7r~ thresholds. 



Table 4.3. Decay widths -Tc/n for both charged and neutral 
decay channel. The corrections for the BKW and DKM ampli- 
tudes are expressed relative to the strong GL result. DKM(uc) 
denotes the variant of the amplitude without the universal pho- 
ton corrections subtracted. 







- 


pGL 


154.5 eV 


222.8 eV 


^pBKW 


(-1.5±1.3) eV 


(-2.5 ± 2.0) eV 




= (-1.0±0.9)% 


= (-1.1 ± 0.9)% 


^jiDKM 


(-2.9±0.7*)eV 


(-7.3 ± 4.0) eV 




= (-1.9±0.5*)% 


= (-3.3 ± 1.8)% 


^jiDKM(uc) 


(-1.5±0.7*)eV 






= (-1.0±0.5*)% 





terminations [55,57,58,68-73].^ Uncertainties due to the 
Ki are on the 1% level throughout. 



The J] Stt" Dalitz plot is displayed graphically for 
our central result, the full one- loop amplitude including 
electromagnetic corrections up to O(^e^), in Fig. 4.4. The 
cusp structures at the tt+tt" thresholds are clearly visible. 
Figure 4.4 also demonstrates that decay spectra with re- 
spect to s, t, u are most sensitive to the cusp effect, not 
with respect to the radial coordinate z. 



4.4 Decay widths, branching ratio, quark mass ratios 

The charged and neutral total decay widths can be calcu- 
lated via 



c/n 



dt \ Ac/nis,t,u)\ 



(4.4) 

where the symmetry factor 5„ = 1/3! for the neutral de- 
cay accounts for the indistinguishable three neutral pions, 
while Sc = 1 for the charged decay. The s-intcgration 



boundaries are given by s™™ = 4M^ 



= 4Mio, and 



c/n 



{Mjj — M^o) , whereas the s-dependent limits for 
the t-integration can be found in Ref. [43]. 

The widths thus calculated are shown in Table 4.3. 
The strong widths calculated at one loop underestimate 
the experimental values significantly, so again, the elec- 
tromagnetic corrections are given as relative corrections to 
the strong result (GL) . For the charged decay the width is 
reduced moderately by about 2%. We also show a variant 
of the DKM amplitude (denoted by DKM(uc)) where the 
universal photon corrections have not been subtracted, i.e. 
the Coulomb pole at s = 4M^ and the kinematical brems- 
strahlung singularity at ,s = (Af^ — Mj^a)"^ arc included. In 
this case, the reduction of the width is only about 1%. 
In the neutral channel, the shifts in the width are com- 
pletely dominated by the corrections to the normalization 
of the amplitude, see Table 4.2. The total corrections up 
to 0{Se'^) are about a factor of 3 larger than those es- 
timated by BKW. The total reduction of the width by 
electromagnetic effects is in good agreement with the one 
found in Ref. [16]. 

From these widths one can determine the branching 
ratio 



(4.5) 



at the different orders in isospin breaking. The results are 

given in Table 4.4. For the BKW corrections, as both 
widths are shifted by almost the same amount, the ra- 
tio is nearly constant compared to the strong result (GL). 
Including all radiative corrections (DKM), r is reduced by 
1.4%, with an uncertainty of about the same size. With- 
out subtraction of the universal soft-photon corrections, 
the effect on r is somewhat larger, about 2.3%. 



* Note that at 0{p^) [13], a can be made consistent with 
experimental results by particular choices of the low-energy 
constants. 
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Table 4.5. Corrections to extractions of the quark mass ratio 
Q from charged and neutral decay channel. For details, see 
main text. 





7} — > n^n n" 




^gBKW 


(+0.24 ±0.22 )% 


(+0.28 ± 0.22)% 


^gDKM 


(+0.48 ±0.12*)% 


(+0.84 ± 0.46)% 


^qDKM(uc) 


(+0.24 ±0.12*)% 





Finally, wc can read off the corrections for the quark 
mass ratio Q from the relation 




(4.6) 



Strictly speaking, the relation (4.6) does not hold for the 
electromagnetic terms of (^(^^(mu = m^)) as present in 
the BKW corrections. However given the smallness of the 
corrections and the size of the uncertainties, the additional 
error in factorizing in the complete amplitude can be 
safely neglected. The resulting shifts in Q are collected 
in Table 4.5. They are to be interpreted in the follow- 
ing way: extracted naively according to (4.6), Q would be 
shifted compared to the purely strong value according to 
Table 4.5, i.e. generally increased. In order to purify a real 
measurement (that includes electromagnetic effects) such 
that the purely strong value for Q is obtained, the opposite 
shift has to be applied. The electromagnetic corrections in 
Table 4.5 are below half a percent for the charged chan- 
nel, while in the neutral channel, a value of Q would be 
reduced by close to 1%. 



5 Conclusions 

We have re-evaluated electromagnetic contributions to 
r] ^ Bit decays, thus calculating the corrections to Suther- 
land's theorem [1,2], extending and updating an earlier 
analysis [14] that neglected terms of C'(e^(m„ — m,d)). 
Terms of this order contain the essential non-analytic 
structures affecting the Dalitz plot distribution, in partic- 
ular real- and virtual-photon corrections in ^ 7r+7r~7r°, 
and the cusps due to virtual tt+tt" intermediate states 
in T] ^ Stt''. We have shown (explicitly at the soft-pion 
point) that these additional terms are not suppressed by 
an additional small isospin-breaking parameter, but only 
by a factor of (m^ — mu)/{md + m„) w 1/3. Terms of 
C'(e^(m„ — rrid)) violate the leading isospin relation be- 
tween charged and neutral decay amplitudes. 

We have calculated the corrections to the various 
Dalitz plot parameters, as well as for the widths and 
branching ratios, and the resulting correction for an ex- 
traction of the quark mass ratio Q. Although the effects 
of 0(e^(m„ — md)) are of the same size as those analyzed 
in Ref. [14] and are not any further suppressed, the total 
electromagnetic corrections remain very small, at the per- 
cent level, throughout. The most significant change occurs 
in the tj ^ 37r° Dalitz plot parameter a, which is slightly 



more strongly affected by the presence of the cusps. The 
latter are, however, not sufficient to explain the sign dis- 
crepancy between chiral predictions and experimental re- 
sults. 

The present study does not replace a dedicated anal- 
ysis of the cusp effect in ?7 ^ Stt" in non-relativistic field 
theory [24,25,74], but rather complements it in also pre- 
dicting the polynomial electromagnetic terms. As such, 
our results ought to serve as correction factors in future 
precision experiments on — > Stt decays. 
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by DFG (SFB/TR 16, "Subnuclear Structure of Matter") is 
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At) ^ Stz amplitudes 



In this appendix, we collect the explicit formulae for the 
full ?7 ^ Stt decay amplitudes. In writing these down, we 
make use of the relation 

Bo6 _ V3{3M^ + M^o) 
' V3(M2-M2„) 16M2„Q2 ' ^ ■ ) 

valid at leading chiral order and to leading order in isospin 

breaking, and use the redundant quantities e, 5, to 
express the various next-to-leading order contributions in 
a more compact form, where possible. For numerical eval- 
uations, only Q is employed, hence e and 5 are always to 
be taken as abbreviations according to (A.l). In a simi- 
lar spirit, we do not perform the expansion of kaon loop 
functions according to (3.10) and (3.11) explicitly, as these 
also tend to make the results more cumbersome. 



A.l Renormalization factors, mixing 



First we write down the Z factors for the wave-function 
renormalization that follow from the full propagators at 
NLO shown in Fig. 3.4 and the renormalization corrections 
defined in (3.18). As the leading-order amplitudes for both 
T] 7r+7r~7r° and r] 37r° are of 0{5), all formulae are 
only needed up to C(e^) in isospin breaking. The Z factors 
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read 



(A.2) 

^ r Ak + Ako + AA^ _ 3 2 
\ 48 2 ^ 



-{M^^+Ml,)L, 



- M^oL5 1 - e^|-{ifi + K2) - 2{2K3 - K^) 
20, 



(A.3) 



7 1,8/ Ak + Ako+2{A^ + A^o) 2 r 



A Q on ^ 

(A.4) 

In order to account for Tr^rj mixing due to the diagrams 
depicted in Fig. 3.2 at next-to-leading order we need the 
quantities £4 and Z^o^ defined in (3.14), 



£4 = 



25oU.o.or , r ^ , MIo{Ak + Ako^2A^) 



V3F2 1 



12(3^7 + is) + 



Ar, - A^o 



{3M^ + M^o){Ak - Ako) 



4(M2 - M^o) j 4x/3F2(M2 - M^^) 



2e^ 



-{2K3-K4)-Kr,-Ke 



3x/3(M2 - M2o) 

x(M^2+M2o)(1-^) 

+ 2{K, + if 10) (m^o - -^(3M2 - M^o)) I 



+ 0(6^ 



(A.5) 



3\/3 V V3 
+ 0(^2), 

that multiply the lowest-order amplitudes 



(A.6) 



A: 



LO 



3s - M2 - 2M2o 



"^r7->r)7r+7r- ~ "^77->777r''7r'' ~ 2"^7r°^37r° ~ ' ^ ') 

where quark mass insertions have been replaced by leading 
meson masses since both Z^o^ and €4 are next-to-leading 
order effects. 



The various renormalization corrections defined in 
Sect. 3.2.2, as well as Aqmo and Ap that we use to ex- 
press strong low-energy constants in terms of observable 
quantities, are given by 



^GMO = ^\ 12{M^ - M^o) [L5 - 6{2Lr + Lg) 
^ \3M^ + M^o){Ak + Ako) 



- Ml, 



2Ml„{2A^ - A^o) - 6M^Ar, 
' Ak-A^ 1 



2e'^{ 3 



M2 - Af 2^ 167r2 



2K3 - Ki 



1 



+ -(if5 + J^6)-6(ifio + ifll) 



2AMl 



Aw = 



F2 

1 f 1 



6L5 



Ak-A^ 
- M2„ 



} + 0(5) , (A.^ 



F3\S 



5A„o - 2Ako - ZAr, 

+ 3(M2-M2o)L5| + 0(e^(^) 

= ^{(3M^' + ^'o)(ZiK + Ako) 
- 2Mlo{4.A^ - A^o) - 2(2M2 - Mlo)A^) 

^ " " ^-{2(M2-M2o)(3L7 + i8) 



(A.9) 



+ 



+ {M^ + M^o) [3(2L6 - U) - 2(L5 - 2^8)] } 



^{2(m2 - Ml,){Ki +K2-Kr- K^) 
-^{M^-3Mlo){2K,-K4) 
+ ^{3M^-5Mlo){K5 + Ke) 
- |(m2 - 2Mlo){K, + if 10)} + 0{S) , (A.fO) 



z = — ^(^TT - A^o) -^^2A^ + Ak 



AMI 



+ 4 



3(M2 + Mlo)U + 2MloL^ } 

- e^^2Ml„ [2if3 -Ki- 4(ifio + ifi 

-Q{M'^+Mlo)K^ + 2,A^]+0{5) , (A.ff) 



1 

8^ 



. 1 \ Ak + Ako . ^3M^-^ArO-A, 
= f^]^ 2 o.,.9 , .^9 



3M2 + Ml, 



2A{M^-Mlo) 



3M2 + Ml, 



M^{Lr, - 2L8) 



{M^ - Mlo){SLr + Ls) 
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4e" 



3M2 + Ml 



+ i(3M2 - M^o){Kg + K,o)} + 0{6) , 



(A.12) 



\^o}+Oie^S). (A.13) 



A. 2 T] —>■ TT+TT tt" decay amplitude 



All loop contributions in the following arc given up to 
the order in isospin breaking considered in this work, i.e. 
0{Se'^). Summing up the tadpole diagrams depicted in 
Fig. 3.5a) for the charged decay channel yields 



itad 



12v/3F4 

+ (8s + - 15M^o)A^o + (M^ + M^o)Ar, 
+ 4(3s - - 5M^o + 4ziM2)z\^| . (A.14) 



Next we give the contributions from the s, t, u channel 
loop diagrams Fig. 3.5b) in terms of the loop functions 
defined in Appendix B. The superscripts s, t, u of the 
amplitudes specify the channel, the subscripts for the t 
channel diagrams the virtual mesons in the intermediate 
state. The s channel contributions are given by 



48V3F, 



^ { 3s(3s - 3M^ - M^„ ) ( Jkk (s) - Ji^o^o (s)) 



+ 



+ 2(5s - 3M2 - MI^){Ak - Ako)] 
^^{(3s - 3M2 - MI,)Jkk{s) + 2Ak\ 

-3{s-M^o)J„o„o{s)-2A^o 

(3s - 4M2o){3sJ^o;^o(s) - SM^o J^^o(s) 

+ 3(s + AAM^){Jkk{s) - 4J,„(s))} 
+ 2(5s - 4M^o)Ako - SMlo{Ar, + A^o) 
+ 2(5s - AMlo + 12AMI){Ak - 4zl^)] } . 



1 

+ 12 



The individual t channel contributions can be written as 

e r 1 



18^4 I 15^2 



-( 



5t+6s-2M^-2M;^o 11 + 3 



{t - &Mlo){t + 2s - M2 - 3Mlo) 



+ 3 



2t^ + t{s - 2M2 - 7M2o) - M2o(4s - 5M^ - TM^o) 

- + 2t{s - M2 - SM^o) - 5M^o - M^o)^ 



AMI 
t 



lt + s-M^^-lM:.-(s^^^t_^\AM? 
X (A^ + A^o) 



3{M^-M^o){t-3M^o) — 



A^-A^ 



(A.16) 



AM'^ ( 

AW = -^71^ I - - Ml.) J Ko Kit) + Ak 
( 3Ml + M%\, J 

+ Ako -[2 ^—^] {^K - Ako) I 

~2^{4^[^*+^'"^^'"^^^°^^^*"^^^'~^^'''^ 
- ^2(5t + 6s - IIM^ - 13M^o) 



_ JM^2-M^,)(3M^2 + M^2e) \ 



4.t^-t{s + 7M^ + Mlo) + {s + M^-Mlo) {3M^ + M^o ) 



- I 7i - 2s - - 3M^o 



(M2-M2o)(3M2+M2o) 



4t - 2s - 5M^ + 



AMf 



jRORit) 



( Ml-Ml.\ 
-2(2+^^-^ — ^jAMl 



iAK + AKo)\,{A.17) 



+ 



"'^ 6Fi I 3 

71 

3 



, M\ ( Ml-M\\ „ 
t-Ml-^-\2^ " ^ ^ \amI 



The total t channel contribution is therefore 
Al''{s,t) = Alo^ + A*KOK + -Al,^ ■ 



(A.18) 



(A.19) 



(A.15) 



The u channel graphs are given by crossing symmetry as 
Ai-{s,u). 
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The triangle diagram shown in Fig. 3.6b) gives rise to 
a contribution of the form 



{3M^-nM^o)AK - 2(3M2+M2o)Zi^} 



2{s-2Ml)G{s) + J^^{s) 



- 6^(35 -4M2o) 



+ 2 



1 



Ml 87r2 



M2 47r2 

(A.20) 



15M2 + Ml„ 



I mi 1 
■ in ■ 



and the vertex correction diagram depicted in Fig. 3.6a) 

leads to 



647r2 Ml 4Ml 3 
(^M^-M^y2Ks-Ki)-M^{Kg-5Kio-QKn) 



_ eM^-4M2) r^_^| +e2(3M,^ + M^o) 

2F2 \M2 127r2/- 



^ + jK-M^o)(2i^3-i^4) 



The result of the various replacements in A^J^'^ for the 
charged decay is given by the total tree amplitude 



A. 



tree '^K - Ako 



(3s - 4M2 )2M2 ^ „ 



- Ml, 



6\/3F4 

e2(3.s^4M2o)2M^o / 26 X 

9V3i^2(M2-M2o) V ^) 

3{2K3 - Ki) - 2{K5 + Ke) + 2{Ks + K^) 
1 



Q22^/3F2^2 



(3s-4M^o) 



3M2 + M2o 



, (3M2 + M2o)M2o+6M2z\M2 

+ 3(M2-M2„) 

3M2 + M2 



M^-M% 



X {(11M2 - 3M2„)^ + (5M2 + TM^o) 

- 3(3Af2 _ M2o)zi^ - {M^ + M^o){A„ + A^o)} 

QAMlMl r / 3e2i^2 
+ ... \4i + 4t^]{Ak-A^) 



A 



M^-Ml, yV ' AMI 
+ 2Ako - 5A„o + 3Z\,, 



-(3M2 + M2o) 

L. 



Ml 



-^GMO + 



{M^ + Ml,)AMl 
M^-Ml, ■ 



^[{s-M^^){s-2AMl)+2tu-{^s+2AMl)Ml,} 



1 /U 



\ M'^,{Ak + Ako) + {M^-AMIo){A^-A^.) 
1 



- Ml, 



4Mloi3M^ - Mlo)Ar, 



Mt + 15M4 



^[Ak + Ako- A, 



A^o A. 



2 V " " 2 
?,e^F^{MdloAK - {3Mlj + 5MIo)A^} 
-3 



Mlo 



- M^Ke + 3{M^+Mlo)iKw + i^ii) 

(A.22) 



A.3 ri — > Stt" decay amplitude 



For the neutral decay channel, the tadpole diagrams de- 
picted in Fig. 3.5a) sum up to 



.tad _ M^{Ak-Ako) e f a , w 

+ 2(m2 - 15Mlo)A^ + 3(3M2 - 5Mlo)A^o 

+ 2(4M2 + 3M2o)(Z\x + Ako)^ . (A.23) 



The rescattcring diagrams displayed in Fig. 3.5b) give 
rise to the following s channel contributions: 



1— {2(5s - 3M2 - mI„){Ak - Ako] 



48\/3F4 



+ 3s(3s - 3M2 - MIo){Jkk{s) - Jkoko{s)) } 
+ -^[{M^-Mlo)Mlo [Jnnis) - 3J,o,o(s) - 2^0(5)] 



12s2 - 2s(9M2 + 5Mlo) + (3M2 + Mlof 
X {Jkk{s) + Jkok°{s)) 



+ -(10s - - ^MIo){Ak + Ako) 

- 2(3s - AMlo){s - Mlo)J„„{s) - ^(5s - 7Mlo)A„} . 

(A.24) 



- ZiM2{-(M2 + m2o)(2Z\^o - 5A^o + 3 An) The crossed channels are given by ^^*"(t) + ^^*"(w). 
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Replacing the strong LECs in A^^^ for the neutral 
decay as explained in the main text yields the following: 



{3(2/^3 - K^) - 2(if5 + Kg) + 2{Kg + i^io)} 
1 f 1 



(^Ak + Ako-A^ + ^-^A^) -^{M'^-Mlo)A^ 

:-A„o-Ar,)^ 



2,Ml-bMl, 



^{Ak-Ako)-7M^o{2A., 

3(9M4 - 2M^M'^o + M^o)iA^ - A^o) 
3 



^GMO , 



AM'z 



- Ml, 



+ 3e 



1 



M2 - M2 167r2 



2(ifio+i^ii) 



{M' - M^o)(3M^ + Mlo)[Ki + K2) 



(A.25) 



- (3M2 - 5M2o)(5M2 + M^o) (i^3 - ^) 
+ M2(3M2 + Ml,){Kr, + Ke) 

- 3M2(Af2 - Mlo)iK9 + if 10) 
B Loop functions 



In this section we give explicit expressions for the relevant 
loop integrals, performed in d dimensions, with a renor- 
malization scale fi. 



B.l Meson loop functions 

We define the basic meson loop functions 



_ 1 /■ d'^fc 1 



i J {2TT)d Ml - fc2 ' 

1 [ cfik 1 



(B.l) 



)<^[M2-fc2][M2-(fc-g)2] 



with s = q . The tadpole function Aa is given by 



Aa = Mi { 2X 



2J.. , _^^Ml 

167r2 /i2 



(B.2) 



with the constant A containing the divergent part in (rf— 4) 

.d-4 ( I ^ 



x = 



167r2 Irf - 4 2 



(ln47r + r(l) + l) 



(B.3) 



where -/"(I) = 7b ~ 0.5772 is the Eulcr Maschcroni 
constant. The loop function Jab{s) can be split up accord- 
ing to 



Jab{s) = Jab{s) - Jab{0) 
Jab{0) = -2A - 2kab , 



(B.4) 



with the scale-dependent constant kab and the function 

Jab{s), 



kab 



1 



Jab{s) = - 



327r2Z\„fc 

1 f 



167r2 



^ab _ ^\ ^Mg ^ 
S Aab) ^ Mb 



+ 



2s 



In ' + - 2ni 0{s - {Ma + Mkf) 



+ - 



S - Sab - V 
\/-v'^ \ , Aab + S Aab-S' 

arctan — , — arctan — , 



(B.5) 



where we have used the abbreviations Sab = M^ + M^ 
and Aab = M^ — M^, the usual Heaviside function 6{x), 
and v"^ =\{s,Ml,M^). 

B.2 Triangle loop function with virtual photon 

For the triangle diagram depicted in Fig. 3.6b) we need 
the integral involving two charged-pion propagators and 
one photon propagator, 



G{s) = 

1 r d'^k 



(B.6) 



I J {2n)d [Ml -(k + qafml -{k- gb)2][fc2 _ ^2] ' 

where s = {qa + qb)^ and we only consider the on-shell 
case q^i^ = Ml. G{s) can be written, in the relevant kine- 

matical region s > AMI above threshold, as 



G{s) = 



1 



327r2scr 
+ 41n 

- 21n 



4 Li 

l-a 
l + a 

ah 

777.2 













' 2o \ 










— a 




h m 




+ (J 



ln2 



1 - cr 



+ OK), (B.7) 



with the dilogarithm or Spence's function 

L.(.)=y^ ^^dt. 



(B.8) 



In this representation the real part of the Coulomb pole, 
i.e. the kinematical divergence at threshold, resides solely 
in the term proportional to — 27r2. 
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C Numerical parameters 

For the meson masses we use M^o = 134.98 MeV, M^r = 
139.57 MeV, Mko = 497.61 McV. Mk = 493.68 McV and 
Mn = 547.85 MeV [43]. For the meson decay constants, 
we employ = 92.2 MeV and Fk = 1.193 [43], the 
latter value relying on Standard Model electroweak cou- 
plings [75]. The electric charge and the quark mass dou- 
ble ratio are evaluated using a = c^/Att — 1/137.036 and 
Q = Qd = 24.2. The photon cutoff energy Ecut is set to 
a typical detector resolution of 10 MeV. 

For the only strong low-energy constant not expressed 
in terms of physical obscrvables, wc use the value L3 = 
—3.5 X 10~^ from Ref. [76]. As we concentrate on the elec- 
tromagnetic corrections in this article, we do not consider 
uncertainties for Z/3. For the electromagnetic low-energy 
constants , wc rely on the estimates in Rcfs. [77 79] 
(compare Refs. [80-82] for alternatives) using the Feyn- 
man gauge and given at a scale /x = Mp = 0.77 GeV: 
-Kl = Kl = 2.7 X 10-3, 4i^5 = 2KI = Kl = 2.8 x lO'^, 
Klj =Kl = 0, K{q = 4.0 X 10-3^ and K^^ = 1.3 x IQ-^. 
Since no number for Kg is offered, wc go back to an earlier 
evaluation in Ref. [80] and use Kg — 0. The uncertainties 
in the K^ are difficult to assess; they are the dominant 
sources of uncertainties for the electromagnetic correc- 
tions. We adopt the following procedure. Since the val- 
ues of the renormalized LECs at two different scales are 
related by 



i^r(M2) = if[(Mi) 



167r2 



In 



Ml 

M2 



(C.l) 



it appears natural to use correlated errors due to a varia- 
tion of scale according to 



KT ^k:± 



Si 
167r2 



(C.2) 



This procedure has the advantage that the resulting er- 
ror estimate is invariant under a redefinition of the La- 
grangian (2.9), in contrast to a more naive uncorrelated 
variation of all KJ according to KJ ± l/167r^. However, 
in some cases (most notably for the normalization of the 
DKM amplitude for 77 tt+tt-tt^), an accidental cancel- 
lation between the various Ei in (C.2) leads to an unreal- 
istically small error. In those cases, we have replaced the 
Si by 1 in (C.2) and marked the corresponding errors, ob- 
tained in a non-standard way, by an asterisk in Tables 4.1, 
4.3, and 4.5. 
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